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In this paper, we calculate the decay rates of D + — > D°e 
Dg — > D + e~e + and Bg — > B°e~e + semileptonic decay processes, in which only the light quarks 
decay, while the heavy flavors remain unchanged. The branching ratios of these decay processes 
are calculated with the flavor SU(3) symmetry. The uncertainties are estimated by considering the 
SU(3) breaking effect. We find that the decay rates are very tiny in the framework of the Standard 
Model. We also estimate the sensitivities of the measurements of these rare decays at the future 
experiments, such as BESTII, super-B and LHC-fe. 

PACS numbers: 13.30.Ce, 13.20.Fc, 13.20.He, 12.39.Hg 



The properties of hadrons containing a single heavy 
quark Q (rug ^> Aqcd) along with light degrees of free- 
dom are constrained by symmetries which are not appar- 
ent in QCD [lj. These symmetries are manifest in the 
heavy quark effective theory (HQET) where the heavy 
quark acts in the hadron's rest frame like a spatially static 
triplet source of color electric field [2], |3( . In the HQET 
the heavy quark's coupling to the gluon degrees of free- 
dom are independent of its mass and described by a Wil- 
son line In the rare decay processes of D + 
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"T" (3 1/ 



D°e+v, 
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B e e , the heavy quark flavors (c or b) remain 
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unchanged, and the weak decays are managed by the 
light quark sectors. In the limit of the flavor SU(3) sym- 
metry of the light quarks, the matrix elements of the 
weak current can be constrained, the uncertainty can be 
estimated. 

In this work, we study the rare decay processes D + — ► 



and B° s 
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v , B v s -> B+e-V, -> D+e-e+ 
Applying the SU(3) symmetry for 



the light quarks, the form factors describing the strong 
interaction in these decays can be obtained. The uncer- 
tainties can be estimated by considering SU(3) breaking 
effect. 

For the semileptonic decays D + — ► D°e + v and D$ — > 
D°e + i>, the decay amplitude is 



element in eq.([T]) can be decomposed as 

(D°(p 2 )\q 1 ^(l- l5 )q 2 \D+ ) (p 1 )) 
= .f+(q 2 Kpi+P2), + f-(q 2 )(Pi-P2)„ (2) 



where f±(q 2 ) are the form factors including all the dy- 
namics of strong interaction. Considering the decom- 
position of the hadronic matrix element in eq.([2]), and 
neglecting the lepton mass, we get the amplitude as 



A = ^V l3 f+(q 2 )( Pl +p 2 )MW(l ~ TBMfa). (3) 



Then the decay width can be obtained as 
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where G F = 1.16639 x 10 -5 GeV is the Fermi constant, 
and Vij is the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix element, the functions u(k±) and v(k 2 ) are Dirac 
spinors, which describe the productions of the neutrino 
with momentum k\ and the anti-lepton with momentum 
k 2 , respectively. 

According to its Lorentz structure, the hadronic matrix 



where mi and m 2 are the masses of the initial and final 
heavy hadron involved in these rare decays. 

Next we consider the flavor-changing neutral current 
(FCNC) processes D+ -> D+e~e + and B° s -> B n e~e + . 
They can be described by the AS = 1 effective Hamilto- 
nian in the quark level at scales n < m c [||: 
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where the operators Qi(n)'s are defined as 
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Here the indexes a and /3 are color numbers, the summa- 
tion ^2 runs over all the quark flavors which are active 
at the scale /z, and the V ± A denotes 7 M (1 ± 75). 

At n = 1 GeV with A$ =215 MeV and NDR scheme, 
the Wilson coefficients Zi and Y^s are calculated to be 
5] Z x = -0.409, Z 2 = 1.212, Z 3 = 0.008, Z 4 = -0.022, 
Z 5 = 0.006, Zq = -0.022, Z 7V = -0M5a QE D, and 
Fx = Y 2 = 0, F 3 = 0.025, F 4 = -0.048, Y 5 = 0.005, 
F 6 = -0.078, Y 7V = 0.747a QE D, Y 7A = -0.700a QED , 
where ckqed is the fine structure constant. 




FIG. 1: Diagrams for the rare decays of D$ — > D + e + e and 
B° s B°e + e~ up to one-loop level, where the cross circles 
denote the operator insertion, (a): the tree level diagram; (b) 
and (c): the one- loop level diagrams with two different ways 
of operator insertions. 



We calculate the FCNC processes D\\ 



and 
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B°e + e up to one-loop level based on the effective 
Hamiltonian. The diagrams we considered are depicted 



in FigJT] The amplitude for the FCNC processes is cal- 
culated to be 

A = [a 1 u{k 1 ) 1 tl v{k 2 ) + a 2 M(fc 1 ) 7 ^7 5 w(/c 2 )] 

x(pi+p 2 W+(g 2 ), (7) 
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where iV c = 3 is the color number, Q u , Q d , Q s and Q e 
are the charge of the relevant quarks and electrons, m s 
is the mass of the strange quark, and the function G is 
defined as 



G(m 2 ,g 2 ,tt) = - I dx4x(l-x)x 
Jo 
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which is originated from the loop calculation. The pa- 
rameter t is defined to be 



v* s v ud 
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With the amplitude in eq.0, the decay width can be 
obtained as 
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In the numerical calculations we use the following values 
for the Wolfenstein parameters 0: 

A = 0.2272; A = 0.818; p = 0.221; fj = 0.340. (13) 

In the limit of SU(3) symmetry for the light quarks 
u, d and s, the masses of the light quarks are the same, 
the masses of the mesons with the same heavy quark but 
different light flavors are also the same. Therefore the 
energy release in the rare decay process is zero in the 
limit of the SU(3) symmetry. The velocity of the heavy 
mesons in the initial and final states is unchanged. The 
wave functions of the light quarks are the same after the 
rare decay occurs. The heavy quark does not feel any 
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change during the light flavor transition process. Thus 
with the flavor SU(3) symmetry holding, the form factor 
f+(q 2 ) can be normalized to unity at the point q 2 = 0, 



/+(0) = 1. 



(14) 



Note that there is only one kinematic value for q 2 with 
the SU(3) symmetry holding, q 2 = 0. However, with 
SU(3) flavor symmetry broken, i.e., the mass of s quark 
m s being larger than that of u and d quarks, the light 
quark is boosted after the light flavor transition occurs. 
The possibility for the light and heavy quarks still bound 
together to form a heavy meson decreases, this leads to 
a form-factor suppression. Thus the value of the form 
factors must be smaller than unity. The larger the s 
quark mass, the larger the deviation. 

With the SU(3) symmetry broken, the allowed value of 
q 2 is extended from q 2 = to a range of < q 2 < <i^ ax . 
At the point q 2 = 9m aa ., all the energy and momentum 
released in the decay are carried away by the lepton pair, 
the velocity of the final meson is not changed (this is the 
point of zero recoil) , only the mass of the light quark 
changed, therefore the deviation from the limit of SU(3) 
symmetry is the smallest at the point of zero recoil. For 
the point q 2 — 0, there is further deviation caused by 
the kinematic change of the final meson. Thus one can 
write the the form factor at q 2 = q 2 ^^ as an expansion 
in terms of a small SU(3) breaking parameter. Including 
the SU(3) breaking effect, eq. (fT4|) is extended to the point 
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where the parameter \su(3) describes the correction due 
to SU(3) breaking effect at q 2 — <?^ a2 .. Note that the 
value of Xsu(3) must be negative, because the SU(3) 
breaking effect leads to a form-factor suppression. In 
addition, if the expansion including the SU(3) breaking 
effect is performed at the point q 2 = instead of the 
point of zero recoil, further correction due to the velocity 
change of the final meson should be considered. 

The larger the mass of the s quark (here we neglect the 
masses of u and d quarks), the larger the SU(3) breaking 
effects. An appropriate SU(3) breaking parameter can 
be taken to be m s /A, where A ~ lGeV is a hadronic 
scale. According to the Ademollo-Gatto theorem 0, [1| , 
the parameter \su(3) should be of the second order of 
the SU(3) breaking parameter, i.e., 



A SC / (3 ) ~ 0((m s /A) 2 ). 



(16) 



In this work the correction to the form factor due to the 
SU(3) breaking effect is treated as an uncertainty, the 
form factor is varied in the range 



l-(m s /A) 2 < f+(q 2 max ) < 1 



(17) 



The momentum transfer squared q 2 should be in the 
range range < q 2 < (mHi — rnH 2 ) 2 , where ran 1 and 



mH 2 are the masses of the initial and final mesons con- 
taining the heavy quark, respectively. The q 2 dependence 
of the form factor is typically governed by the nearest res- 
onance, which can be parameterized as a pole-dominance 
form 



U(q 2 ) 
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1 - q 2 /m 2 H , 



(18) 



where H* is nearest pole resonance which contributes 
to the from factor. H* should be K* in s — > d, or u 
transitions. With eqs. p"?")) and (|18p . we can obtain that 
the form factor at q 2 — should be varied in the range 

(\-{m s /h) 2 )(l-q 2 max /m 2 H ,) < /+(0) < (\-q 2 nax /m 2 H ,). 

(19) 

Taking the mass of strange quark m s = 100 MeV, by 
considering the SU(3) breaking effect, we get the follow- 
ing branching fractions: 
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= (6.44 ±0.07) x 10~ 17 



(20) 



where the uncertainties are caused by SU(3) breaking 
effect. For D+ — > D°e+i', the correction to the form 
factor due to isospin breaking effect is neglected in 
the above presentation, the numerical correction due to 
isospin breaking should be far less than the effects in the 
strangeness- violating processes. Here we take m c = 1.5 
GeV, raj, = 4.8 GeV in the numerical calculations. 

The mass of the strange quark is very close to the mass 
differences of the initial and final mesons in the decay pro- 
cesses involving s quark, D+, — > D°e+v, D+, — > D+e+e~ , 
B° s -> B+e-v, B% -> S°e+e", which is about 100 
MeV. It is interesting to study the role of m s in these 
decay processes by slightly varying the mass of strange 
quark. The results are given in Table [U The decay rates 
are only slightly changed with m s varying in the range 
100 MeV <m s < 140 MeV. 

Within the SM framework, we find that the branch- 
ing fractions for these rare decays are tiny. However, 
in the coming experiments at BES-III 0, LHC-6 
and super-B factory [l^, it is interesting to search for 
these semileptonic decays. Especially, the decays of 



d: 



~v and Bg 



B+e~v may be reached at 
the super- B factory. 

Since the electron is very soft, one cannot recon- 
struct both the electron and neutrino in the experiment 
near the charm meson threshold at e+e~ colliders. At 
BES-III, to search for the decay D+ — * D°e+v on the 
V>(3770) peak, the charged D mesons are produced in 
pairs, e+e~ -> -0(3770) -> D+D~. Thus, the following 
six tag modes, D~ — > K+ir~7r~, K+7r~ir~Tr° , Ksir~, 
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TABLE I: The decay rates of the rare processes involving strange quark. The second uncertainty is estimated from SU(3) 
breaking effect due to the mass of the strange quark. 
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m s = 120MeV 


m s = 140MeV 




(2.97 ±0.03) x 10~ 8 


(2.96 ±0.05) x 10~ s 


(2.94 ±0.06) x 10~ s 


Be; -> B+e'v 


(4.46 ± 0.05) x 10~ s 


(4.43 ±0.07) x f0~ s 


(4.41 ±0.09) x 10" 8 


D+ -> D+e + e~ 


(3.56 ±0.04) x 10" 1 ' 


(3.64 ±0.06) x 10" 1 ' 


(3.63 ±0.07) x 10" 1 ' 


B» s -» S u e+e- 


(6.44 ±0.07) x 10" 1 ' 


(6.40 ±0.10) x 10 _1/ 


(6.37 ±0.13) x 10" 1 ' 



KsTT~Tr~Tr + , Kstt~tt° and K + K~tt~ , can be used to 
fully reconstruct one of the charged D mesons. The 
summed branching fractions of the six tag modes are 
about 28% of all the charged D decays [6J. The tag ef- 
ficiency for the charged D mesons are about 20%, which 
means that 20% of all the D + D~ pairs can be tagged [Io| . 
For this case, in order to detect the decay D + — ► D°e + i', 
one can reconstruct the neutral D meson decay by using 
46% of all of the D° decays modes in the tagged charged 
D sample [13]. If we see any event of the production of 
the neutral D mesons against the charged D mesons, it 
indicates the observation of the rare semileptonic decay. 
For the decay Dg — > D a e + v, the same method can be 
applied by using the data collected at the center of mass 
E CM = 4170 MeV. With 20 fb" 1 data on the ^(3770) 
peak, the sensitivity of the measurement of D + — > D°e + v 
can be 10~ 6 at the BES-III experiment, while, for the 
measurement of Dg — > D°e + i', it can reach 10 -5 level 
with 20 fb" 1 data running at E C m = 4170 MeV. These 
estimations are listed in table [III 

These semi-leptonic decays can also be searched in the 
B and super-i? factories by using data at T(AS) peak. 
One can reconstruct the following decay chain to search 
for the rare decays D + — > D°e + v: 
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(21) 



where the D* + is boosted, and both tt° and electron could 
have momentum with a few hundred MeV, which can be 
detected and reconstructed in the detector. Since the 
missing neutrino has very low momentum, one can par- 
tially reconstruct the decay of D + , and looking at the 
mass difference Am — m((D° 'e + )w® f t ) — m(D°e + ). The 
signal events should peak around the mass difference of 
m D ,+ — m D + = 140 MeV on the Am distribution. In 
the mass difference, the uncertainty of the reconstruc- 
tion of (D°e + ) can cancel. The resolution on the Am 
will be dominated by the detection of the soft ir°. This 
is a powerful variable to separate background from the 



signal events. For the decay D~g 
the reaction: 



v. one can use 
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(22) 



to extract the rare decay signal by looking at the mass 
difference Am = m{{D°e + )~f soft ) - m((D°e+). With 1 
ab _1 and 50 ab _1 luminosity at B factories and super- 
B, the sensitivities could be 10 -8 and 10~ 10 respectively. 
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In table [TTJ the sensitiveities of the measurements of the 
rare D + and Dg decays are summarized at different ex- 
periments. 

At super-B factory, the data taken at T(bS) can be 
used to search for the rare decays Bg — > B + e~D and 
; + . The cross section of the T(5S) produc- 
collisions is a(e+e- -> T(5S)) = 0.301 ± 
0.002 ± 0.039 nb 0. Unlike the T(45) state, T{5S) is 
heavy enough to decay into several B meson states, in 
which the vector- vector (B* S B* S ) and vector-pseudoscalar 
(BgBs + BgBg) combinations are dominant 1141. About 
30% of the T(5S) decays into B s final states [1^. With 
30 ab _1 data at T(5S) peak at super-B, the sensitivity 
of the measurements of the rare Bs decays can be 10 -9 
by assuming 30% efficiency for the Bs reconstruction. 



TABLE II: Experimental sensitivities at BES-III, B factory, 
Super- B and LHC-b for the rare decays. We assume the inte- 
grated luminosities are 20 fb" 1 (BES-III at ^(3770) peak and 
4170 MeV), 1 atT 1 and 50 atT 1 at B factory and Super-B 
(at Y(4S) peak), 10 fb" 1 at LHC-fe, respectively. 



Decays 
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B factory 
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LHC-& 




(xlO" 6 ) 


(xlO -8 ) 


(xlO" 10 ) 


(xlO" 9 ) 


D + _ D " e + V 


1.0 


1.1 


2.3 


3.8 


D + _» D° e + V 


5.0 


1.1 


2.3 


3.8 


D+ ->• D+e+e' 


11.5 


2.0 


4.6 


5.0 



All of these rare decays can be studied at the LHC-6, 
in which production cross section for bb is 500 About 
5 x 10 12 bb pairs will be obtained in 10 fb -1 integrated lu- 
minosity in the first five years running at the LHC-6 11 1 . 
Accordingly, about 40% (40% or 10%) of the bb pairs is 
predicted to form B° meson (B + or B s meson). Thus, 
the sensitivity of the measurements of the rare Bs de- 
cays is estimated to be 10 -10 at the LHC-6 with 10 fb -1 
luminosity. Both the D* + and D* s + can be reconstructed 
in the decays of B and B s mesons, the decay chains in 
Eqs. |2"T1 and |2"21 can also be used to extract the rare decay 
signals of D + and D~g mesons. The estimated sensitiv- 
ities for the rare decays of the charm mesons are about 
10" 9 at the LHC-6. However, one has to be careful that 
there are many unexpected backgrounds at the LHC-6 
experiment. The estimated sensitivity should be more 
conservative. 

In summary, for the first time, we calculated the de- 
cay rates of the rare D + , Dg and Bs decays, in which 
only the light quarks decay weakly, while the heavy fla- 
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vors remain unchanged. Applying the SU(3) flavor sym- 
metry, the form factors describing the strong interac- 
tion in these decays can be obtained. Considering SU(3) 
symmetry breaking, the uncertainty for the form factors 
can be estimated. Therefore, these rare decays can be 
predicted with the uncertainty estimated by considering 
SU(3) symmetry breaking. We also estimated the sen- 
sitivities of the measurements of these rare decays at 
the future experiments, such as BES-III, super-JS and 
LHC-fe. Especially, the decays of Dg — ► £>°e + f and 
Bg — > B + e~v may be reached at the super-B factory. 
Observations of these decays will be used to test the 
SM predictions for the rare simleptonic decays. Further 
more, any indication of deviation from the SM prediction 
may shed light on the searches of New Physics. 
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